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ABSTRACT 

< 

^ ' We measured polarized dust emission at 350/im towards the high-mass star 

^ ■ forming massive dense clump IRAS 20126+4104 using the SHARC II Polarime- 

ter, SHARP, at the Caltech Submillimeter Observatory. Most of the observed 

2 ■ magnetic field vectors agree well with magnetic field vectors obtained from a nu- 

.^.^ ' 

c/2 ■ merical simulation for the case when the global magnetic field lines are inclined 

c3 ■ 
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with respect to the rotation axis of the dense clump. The results of the numeri- 
cal simulation show that rotation plays an important role on the evolution of the 
massive dense clump and its magnetic field. The direction of the cold CO 1—0 
bipolar outflow is parallel to the observed magnetic field within the dense clump 
as well as the global magnetic field, as inferred from optical polarimetry data, 
indicating that the magnetic field also plays a critical role in an early stage of 
massive star formation. The large-scale Keplerian disk of the massive (proto)star 
rotates in almost opposite sense to the clump's envelope. The observed magnetic 
field morphology and the counter-rotating feature of the massive dense clump 
system provide hints to constrain the role of magnetic fields in the process of 
high mass star formation. 

Subject headings: polarization — techniques: polarimetric — stars: formation — 
ISM: clouds — ISM: magnetic fields — submillimeter: ISM 



1. Introduction 



Massive stars deposit energy into the interstellar medium, hence they may play a key 
role in regulating star formation in galaxies. Nevertheless, massive star formation is still 
poorly understood. 

There are many difficulties in studying the formation of massive stars. Massive stars are 
rare and evolve quickly. The youngest phases of massive stars is poorly understood because 
they are deeply embedded in dense clumps. Often times massive stars form in cluster regions, 
which makes it difficult to disentangle the activities of each young massive star. 

The early B type massive (proto)star, IRAS 20126-1-4104 (~ 7 M©), is a well studied 
unique target that allows us to observe the early stages of massive star formation in a simple 
configuration. The natal massive dense clump, which is a site of massive star formation, is an 
isolated rotating (2 k nis~^pc~^') clump (IShinnaga et al.ll2008l ) with a mass of ~ 200 Mq and 
temperature of 40 K (Shepherd et al.ll2000f). It has a larg e scale bipolar outflow with a size 



scale of 0.5 pc (IShepherd et al. 



2000 



Shinnaga et al.ll2008[ ). The direction of this flow is very 



different from the dir ection of the smaller scale jet emanating from the massive (proto)star 
(ICesaroni et al.lll999l ). indicating that the jet might be precessing. The radial column density 
profile of the massive dense clump shows a shallow slope (r""'^) in an inner region within a 
radius of ~ 0.13 pc, while the outer region of radius > 0.13 pc has a steeper slope (~ r~^-^), 
indicating that the inner region may be experiencing infall, while the infalling wave has 
not yet reached the outer region (IShinnaga et al.l |2008| ) . Another interesting feature of the 
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massive dense clump is tha t the Keplerian disk, known to be associated with the (proto)star 
( ICesaroni et al.lll999ll2005l ). rotates almost in an o pposite sense with respect to the rotation 
of the massive dense clump (jShinnaga et al.l |2008| ) . 



Magnetic fields may play critical roles in star formation ( jMcKee fc Ostrikerl 120071 ). We 
here report a study of the magnetic field of the massive dense clump, IRAS 20126+4104. 



Observations and Data Reduction 



We used the SHARC II Polarimeter, SHARP flLi et al.M2008r ). with the 10.4 meter 
Leighton telescope at the Caltech Submillimeter Observatory (CSO) to measure dust polar- 
ization in the r nassive dense clump at 350 /xm during 2007 August. SHARP uses the detector 
in SHARC II JPowell et al.l l2003h which has 384 pixels. Using the optics of SHARP, one 
divides the detector array into two sections, referred to as the H and V subarrays, for two 
orthogonal polarization components. Both polarization components are thus observed simul- 
taneously. Each H and V subarray has 12 x 12 pixels with a field of view of ~ 1' x 1'. The 
beam size was measured to be 9". The data were acquired under good weather conditions. 
The 225 GHz opacity was measured to be from 0.04 to 0.07. Using the peak of the dust 
continuum emission, pointing corrections were made. Pointing errors were estimated to be 
under ~ 2". A reduced- x^ analysis of Q and U data (JDavidson et al.ll201ll ) was performed 
to obtain the final polarization map. The reduced x^ was found to be 1.87 ± 1.21 when 
dividing the data into eight bins. The error bars for the polarization degree and magnetic 
field direction were inflated accordingly. 



3. Results and Discussion 



3.1. Measured Polarization and Magnetic Field 

Figure [T] (a) shows the polarization vectors measured at 350 /im with SHARP. The 
degree of polarization in the envelope region is higher than that in the central region of 
the clump. Such a ter idency, i.e., polar ization h ole, is observed toward s other star f o rming 
regions such as Orion ( 1Schleuninglll998l ) and W3 (ISchleuning et al.ll2000l ). ISchleuningI (119981 ) 
argues that the decrease of polarization degree may be caused by temperature and/or optical 
depth effects. For the case of IRAS 20126+4104 , the temperature structure appears to be 
complicated, as discussed in lShinnaga et al.l ( 120081 ). The fact that vectors outside of the third 
contour from the peak of the dust continuum have larger polarization degrees compared to 
the vectors inside the third contour (Figure 1) suggests that the polarization hole may be 
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caused by optical depth effects. iMatthews et al.l (120091 ) reported the dust polarization of 



this object at 850 yum. The polarization vectors measured with SHARP are located in the 
central r ~ 0.3 pc region of the massive dense clump, while the polarization vectors measured 
at 850/im trace a region outside of r ~ 0.2 pc but within the r ~ 0.5 pc region. Overall, 
the polarization vectors of Matthews et al. (2009) agree well with our 350 ^m polarization 
vectors. 

The magnetic field directions are obtained by flipping the polarization vectors by 90° 
as plotted in Figure [T] (b). Looking at the magnetic field directions along a north-south 
line passing through the center of the massive dense clump, one notices that many of these 
vectors (except for the central ~ 0.1 pc region) roughly follow the north-south direction. 
On the other hand, for the east and west sides of the clump, many of the magnetic field 
directions tend to follow an east-west direction rather than a north-south direction. This 
tendency persists in the regions where 850;um polarization vectors are observed. Detailed 
comparison between the 350/im polarization vectors and the 850/im polarization vectors will 
be described elsewhere. The magnetic field changes its direction inside the infalling region, 
i.e., for radii below 0.1 pc. 



3.2. Comparison with Simulation Results 

To investigate the morphological evolution of magnetic field lines, we calculated the 
evolution of a magnetized cloud using a three dimensional resistive magnetohydrodynamics 
nested grid code with an isothermal equation of state and sink cell treatment (Machida et 
al. 2005; 2011a). Details are described in Kataoka et al. (in preparation). To conduct the 
simulation, we set the parameters so that they are close to the observed parameters of the 
massive dense clump IRAS 20126-1-4104. 

The schematic diagram in Figure [2] shows the initial state, which is a spherical cloud core 
with a critical Bonnor-Ebert (BE) density profile. The BE profile is characterized by two 
parameters, the central number density [ric) and the isothermal temperature (Tiso). We adopt 
Tic = 6 X 10^ cm~^ with the density enhancement factor of f=1.68 (Machida et al. 2011b) 
and riso= 40 K. Then, a uniform magnetic field {B = 1.5 x 10^^ G) parallel to the z axis 
is imposed on the whole computational domain. The reason why we set the magnetic field 
direction to be parallel to the z axis is because the direction of the observed global magnetic 
field is almost north-south (see Section 3.4). In addition, rigid rotation [Qq = 1.1 x 10"^'^ s~^ 
= 0.35 km s~^ pc^^ ) is added to the initial state, in which the rotation axis, defined by 
right-hand rule, is inclined with respect to the magnetic field (i.e., z-axis) at an angle of 60° 
on the y-z plane, as shown in the schematic diagram of Figure [2J The observer views the 
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simulation from a point in the x-y plane. This is shown in Figure 2, where it can be seen 
that the observer's hne of sight is inchned by an angle of 30° with respect to the y axis. 
On the plane of the sky, the rotation axis becomes P. A. ~ —40°, close to the P. A. = —35° 
rotation axis of the large-scale Keplerian disk (radius of 7400 AU (0.037 pc); Cesaroni et 
al. 1999), which is situated at the center of the dense clump. With these assumptions, the 
initial cloud has a radius of 1.2 x 10^ AU (= 0.6 pc) and mass of 78 Mq. The initial cloud has 

tne energy ratios OI t^thermal/ ^gravity ^ U.O, i^rotation/ ^gravity ^ U.U/ and t^rnagnetic/ ^gravity ^ 

0.55, where ^thermah ^rotation, ^magnetic, and ^gravity are thermal, magnetic, rotational and 
gravitational energy, respectively. In the calculation, we assumed the protostar formation 
to occur when the central density exceeds the threshold number density Uthr = 10® cm~^ in 
the region of radius less than accretion radius, lacc, where lacc = 64 AU is adopted as the 
sink radius (for details, see Machida et al. 2011a). With this treatment, we calculated the 
cloud evolution until the gas accretion onto the protostar or circumstellar disk almost halts. 
To directly compare simulation results with our observations, we calculated polarization of 
thermal dust emission using the formulation in Tomisaka (2011), in which we assume that 
the whole region is optically thin and isothermal for simplicity. 

The bottom diagram of Figure [2] shows our simulation result, in which magnetic field 
vectors and column density 7.1 x 10^ yr after the cloud collapse begins (or 8 x 10^ yr after 
the protostar formation) are plotted. The 30° viewing angle (the top diagram of Figure [2]) 
was chosen because it gives the best agreement between the observed and simulated magnetic 
field vectors. In particular, note that both Figure 1 (b) and Figure 2 show an S shape as 
one moves from north to south. At the epoch shown in the figure, the protostar has a mass 
of 7.3M0 and a large fraction of the cloud mass remains in the infalling envelope. Note that 
radiative heating from the central protostar is ignored. 

In a real high mass star forming region, one can expect in general that the rotation axis 
of a dense clump will not be aligned with the magnetic field. In this situation, simulations 
sometimes show the magnetic field morphology in an hourglass shape, but not always. Also, 
the observed shape of the magnetic field changes depending on the observer's viewing angle. 
For the simulation shown in Figure [2l the observed magnetic field vectors are aligned nearly 
in an S shape, but the morphology of the magnetic field vectors takes an hourglass shape if 
one observes the magnetic field along the x axis. 

In order to reproduce the observed magnetic field morphology, the effect of the rotation 
of the dense clump was found to be essential. This may imply that the rather fast rotation 
of the massive dense clump may have determined the magnetic field directions within the 
central regions of the massive dense clump, as well as the axis of the jet /Keplerian disk 
system associated with the central massive (proto)star. In our simulation, we do see the 
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expected toroidal field generated by the rotating disk at the cloud's center. We cannot 
resolve it in Figure 2 because the scale of the toroidal field (~ 10^ AU) is much smaller than 
the cloud scale (~ 10^ AU). Although the dominant components of the observed magnetic 
field presumably come from large scale structure, the innermost vectors of our magnetic field 
map may show a hint of the toroidal disk field. Although the dominant components of the 
observed magnetic field come from large scale structure, we might have detected a hint of 
toroidal field generated by the rotating disk. 

The configuration of the magnetic field lines is related to the rotating disk that has 
a size of ~ 1000 AU. Since the sink radius adopted in this study (64 AU) is much smaller 
than the rotating disk, we can spatially resolve it. However, the resultant disk size may be 
dependent on the sink radius because we cannot resolve the early phase of the disk evolution. 



3.3. The Magnetic Field and Its Relationship with the Clump's Rotation, the 

Bipolar Outflow, and the Jet Directions 



Based on the observed velocity gradient of the narrow line components, the rotation 
axis of the massive dense clump is roughly 140 ± 20° in P. A., where we again use the right- 
hand rule (jShinnaga et al.ll2008l ). On the other hand, the edge-on Keplerian disk of the 
massive (proto)star appears to rotate in almost the opposite sense to the rotation of the 
envelope of the mass i ve dense clump whos e effective radius is 0.5 6 pc (=1.15 x 10^ AU) 
(ICesaroni et al.ll2005t IShinnaga et al.l |2008[ ) . iMachida et al.l (120061 ) discuss the cases when 
the rotation axis of a circumstellar disk is inclined with respect to the magnetic field axis. 
Under the conditions considered in their paper, i.e., a quiescent cloud core that forms low- 
mass stars, the counter rotatio r i betwe en cloud core and disk, as seen in IRAS 20126+4104, 
cannot occur. iMachida et al.l (l2011bl ) discuss the importance of magnetic braking on the 



circumstellar disk formation in a strongly magnetized cloud. The kind of counter rotation 
observed in this object might happen when the angular momentum of a circumstellar disk 
is extremely efficiently transferred through magnetic braking. 

The two diagrams of Figure [3] show the cold bipolar outflow traced with CO J = 1 — 
and the directions of the jet/Keplerian disk plane, overlaid on the magnetic field directions 
and the 350 //m dust continuum map (same as Figure [T](b)). The jet axis has P. A. ~ —60°. 
Comparing the bipolar outflow axis with the directions of the magnetic fleld in the massive 
dense clump, they are nearly parallel. For the warm bipolar outflow traced with CO J = 6—5, 
the direction is very different (see Figure 12 of Shinnaga et al 2008), particularly the blue lobe 
is in a northwest-southeast direction with a smaller size than that of CO 1 — 0. The direction 
of the CO 6 — 5 outflow is close to the direction of the jet emanating from the massive 
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(proto)star. The observed magnetic field directions near the Keplerian disk appear to he 
paraUel to the disk plane, which means at this location the fields are nearly perpendicular 
to the jet axis. 



Ciardi fc Hennebelld ( 120101 ) report numerical simulations of outfiows in collapsing dense 



cores with misaligned rotation and magnetic field axes and find that a larger angle a between 
the rotation axis and the magnetic field direction leads to decreased efficiency of mass ejection 
via outfiows. The reason why the size of the CO J = 1 — bipolar outfiow in IRAS 
20126+4104 is larger than that of the CO J = 6 - 5 outfiow fishinnaga et al. 



that a for the 1 — outfiow is smaller than that of the 6 — 5 outfiow. 



2008h mav be 



Ciardi fc Hennebelle 



( 120101 ) also found that the misalignment of magnetic field and rotation axes leads to jet 
precession. The apparent precession of the jet associated with the massive (proto)star in 
IRAS 20126+4104 thus may be explained by the misalignment of magnetic field and rotation 
axis. An alternative explanation for precession is that the source at the clump center is 
binary, but observational studies carried out to date do not show evidence of binarity. 



3.4. The Magnetic Field in the Massive Dense Clump and the Global 

Magnetic Field 

It is of interest to compare the magnetic field direction in the massive dense clump and 
the global magnetic field surrounding the clump (e.g., Li et al. 2009). The global magnetic 
field direction around this object is estimated using the optical polarization data archive of 
Heiles (2000). The polarization data points that meet the following criteria were selected to 
trace the global magnetic field: (1) within ± 200 pc from the distance of the object (1500 
pc) for the direction along the line of sight, and (2) within 2° (roughly 50 pc in radius) from 
the center of the massive dense clump on the plane of sky. The polarization data points 
and the mean direction of the global magnetic field are plotted in Figure HI Note that four 
polarization data points that are within ~ 0.8° of the center of the nearby supernova remnant 
(SNR) G78. 2+2.1 are omitted in this diagram because their polarization vectors are likely to 
be affec ted by the SNR . The threshold of 0.8° was determined based on the radio continuum 
map of IWendkerl (119841 ). One sees that the global magnetic field is more or less aligned in 
a north-south direction, i.e., in the same direction as the CO 1-0 bipolar outfiow. The 
massive dense clump is elongated in a north-south direction. The mean direction of the 
global magnetic field is P. A. ~ —3°. The fact that the CO 1-0 bipolar outfiow is parallel to 
the global magneti c field suggests that magnetic field had a significant infiuence on the star 
formation process (JMatsumoto fc Tomisakall2004l ) . However, the deeply embedded Keplerian 
disk and the jet show no such alignment with the global magnetic field. They seem to be 



aligned with the rotation axis of the dense clump. 



4. Summary and Future Work 

In order to measure the magnetic field structure of the massive dense clump IRAS 
20126+4104, polarized dust emission was measured at 350yum. The observed magnetic field 
is consistent with the magnetic field vectors of a simulation when magnetic field lines are 
inclined from the rotation axis of the dense clump. As an evidence that the magnetic field 
played a critical role on the formation of the massive (proto)star, the magnetic field directions 
within the massive dense clump are parallel to the CO 1 — bipolar outfiow direction and 
appear to connect to the global magnetic field. Our observational data combined with the 
simulation results indicate that the rotation of the massive dense clump affects the magnetic 
field within the clump as well as the orientation of the Keplerian disk/jet system. The 
counter-rotation feature observed between the envelope of the dense clump and the large 
rotating disk associated with the massive (proto)star might be a result of efficient angular 
momentum transfer via magnetic braking (e.g. Machida et al. 2011b). This issue should be 
further addressed in the future in order to fully understand the role of magnetic field in the 
star formation processes. It is important to measure the magnetic field strength within the 
massive dense clump using Zeeman effect in order to investigate the details of the evolution 
of the massive (proto)star and of the natal dense clump. 
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Fig. 1. — Black bars overlaid on the 350/xni dust continuum emission (grayscale) represent 
(a) the measured polarization E vectors at 350//m and (b) the measured magnetic field 
directions. The length of the black bars in (a) is set to be proportional to the polarization 
degree. The length corresponding to 5 % polarization degree is shown at bot tom right of the 
diagr am. The contours represent the 350/im image obtained with SHARC II fjShinnaga et al. 
20081 ) and are drawn at 2a, 5a, 9a, 27a, 81a and 243a", where 2a corresponds to 200 mJy 
beam~^. Thick bars and thin bars are vectors with signal-to-noise ratio between 2.5 and 6a 
and between 2 and 2.5a, respectively. 
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Fig. 2. — Top: The schematic diagram summarizes the initial state of the simulation. The 
blue sphere represents the spherical cloud core, being penetrated by magnetic field (purple 
vectors) along the z axis. The rotation axis of the cloud core is inclined by 60° from the z 
axis on the y-z plane, as shown by the big red arrow. Note that observers view the cloud 
core from a direction that is inclined by 30° from the y axis on the x-y plane. Bottom: 
The diagram shows the resultant magnetic field vectors (black bars) and the column density 
(grayscale and white contours) calculated from the numerical simulation. White contours 
and grayscale represent the column density of the dense clump. The protostar has a mass 
of 7.3M0 and a large fraction of the cloud mass remains as the infalling envelope. The size 
of the diagram is 1.2 xlO^ by 1.2 xlO^ AU (=0.6 pc). 
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Fig. 3. — L eft: Molecular b i polar outflow lobes traced in the CO J = 1 — line (same as 
Figure 13 of IShinnaga et al.l (|2008[ )) overlaid on the diagram of Figure 1 (b). Right: White 
line and light blue arrows represent the direction of the Keplerian disk {P. A. -^ 53°) and the 
direc tion of the jet {P. A. ~ — 60°) associated with the massive (proto)star (ICesaroni et al. 
19991 ). marked with yellow star, overlaid on the diagram of Figure 1 (b), magnified in the 
central region. 
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Fig. 4. — Optical polarization vectors (black thick lines) are superposed on an optical image 
(DSS 0.5//m, gray scale) centered on IRAS 20126+4104. The yellow star marks the position 
of the massive (proto)star. Purple thin lines represent the mean direction of the global 
magnetic field [P. A. ~ — 3 deg). Blue thick dashed arrows show the direction of the CO 
1 — bipolar outflow {P. A. ~ — 10°). The red arrow and the light blue arrows represent 
the rotation axis of the large Keplerian disk {P. A. ~ — 37°) and the directions of the jet 
{P. A. ~ — 60°) associated with the massive (proto)star, respectively. 



